Shedding light on the nature of ESO 243-49 HLX-1 
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ABSTRACT 

The point-like X-ray source HLX-1, close to the SO galaxy ESO 243-49, is one of the strongest 
intermediate-mass black hole candidates. We discuss the hypothesis that ESO 243-49 is under- 
going a minor merger with a gas- rich disc galaxy. We propose that the counterpart of HLX-1 
coincides with the nucleus of the secondary galaxy. We re-analyze the available photometric HST 
data, and we compare them with the results of N-body /smoothed particle hydrodynamics simula- 
tions. In particular, we derive synthetic surface brightness profiles for the simulated counterpart 
of HLX-1 in six HST filters, ranging from far ultraviolet (FUV) to infrared wavelengths. Such 
synthetic profiles include a contribution from the stellar population associated with the simulated 
disrupted satellite and a contribution from an irradiated disc model. These are in agreement with 
the observed surface brightness profiles of the HLX-1 counterpart, provided that the merger is at 
sufficiently late stage (> 2.5 Gyr since the first pericentre passage). The main difference between 
models and observations is in the FUV band, where the HST image shows a fuzzy and extended 



emission. 



Subject headings: galaxies: 
X-rays: individual: HLX-1. 

Introduction 



interactions - methods: numerical - galaxies: individual: ESO 243-49 



The point-like X-ray source HLX-1, with a 
maximum luminosity ^ 10^^ erg s^^ (Farrell et 
al. 2009), is the brightest known ultraluminous 
X-ray source (ULX) and one of the strongest 
intermediate-mass black hole (IMBH) candidates. 
HLX-1 is located in the outskirts of the SO galaxy 
ESO 243-49 (luminosity distance ~ 96 Mpc), - 
0.8 kpc out of the plane and ^ 3.3 kpc away from 
the nucleus. The nature of the HLX-1 counterpart 
is still debated. Fits to the Hubble Space Telescope 
{HST) data (Farrell et al. 2012) indicate a total 
stellar mass of < 6 x 10^ Mq. Both a young (^ 10 
Myr) and a old 10 Gyr) age are possible. 



It is not easy to explain how a ^ 10* IMBH 
can form in a 10^ star cluster (SC). Re- 
peated mergers between stellar black holes and 
other compact objects and/or stars were proposed 
to produce IMBHs in globular clusters (e.g. Miller 
& Hamilton 2002), while runaway collapse of mas- 
sive stars can occur in young massive SCs (e.g., 
Portegies Zwart & McMillan 2002). An alterna- 
tive scenario predicts that the IMBH is associated 
with the nucleus of a satellite galaxy which is being 
disrupted in a minor merger with ESO 243-49 (So- 
ria et al. 2010, 2012; Webb et al. 2010; Mapelfi et 
al. 2012, hereafter M12). In this case, the IMBH 
would belong to the low-mass tail of the distri- 
bution of super massive black holes (SMBHs), lo- 
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cated at the centre of galaxies. There is increas- 
ing evidence of galaxy nuclei hosting SMBHs with 
mass < 10^ Mq (e.g. Filippenko & Sargent 1989; 
Barth et al. 2004; Greene & Ho 2004 Secrest et al. 
2012, and references therein). The hypothesis that 
ESO 243-49 recently underwent a minor merger is 
consistent with various features, such as the pres- 
ence of prominent dust lanes around its nucleus 
(e.g., Shabala et al. 2012) and the evidence of UV 
emission centred on its bulge, suggesting ongoing 
star formation (Kaviraj et al. 2009). 

In the current paper, we investigate the minor- 
merger scenario, by analyzing a high-resolution 
A^— body/smoothed particle hydrodynamics (SPH) 
simulation, and by comparing the simulation out- 
puts with the HST data. In particular, we com- 
pare the surface brightness profiles (SBPs) esti- 
mated from the simulation with the observations, 
from infrared to ultraviolet (UV) bands. 

2. HST data revisited 

We re-analyzed the HST data, presented by F12 
(GO program 12256). The data were acquired 
with the Advanced Camera for Surveys (ACS) 
Solar Blind Camera (SBC) in the F140LP fiher 
(far UV, FUV), and with the Wide Field Cam- 
era 3 (WFC3) in the F300X (near UV, NUV), 
F390W (C), F555W (V), F775W (I), and F160W 
(H) filters. The FUV, C, I color-combined im- 
age of ESO 243-49 is shown in Fig. [H It shows a 
fuzzy emission in FUV associated with the HLX-l 
counterpart. Two 'tails' departing from the source 
and leading south and west are visible in the im- 
age. One of the tails seems to be connected to 
a galaxy which was identified with a background 
galaxy at z ^ 0.03 by Wiersema et al. (2010). 
This fact may lead to the conclusion that there is 
some connection between the counterpart of HLX- 
1 and the background galaxy. However, such a 
direct connection seems unlikely because the only 
line visible in the optical spectrum of the counter- 
part (interpreted as Ha) indicates a redshift for 
the HLX-l counterpart closer to that of ESO 243- 
49 (z = 0.0224) than to that of the background 
galaxy {z ^ 0.03, see also Soria & Hau 2012). 

The SBPs in 5 different filters, from F140LP 
to F775W, were derived performing photometry 
within circular apertures of increasing radii (up 
to 1 arcsec), and computing the background in 




Fig. 1. — Color-composite image of ESO 243- 
49. Red is F775W, green is F390W, blue is 
F140LP. The arrow indicates the HLX-l counter- 
part. North is up and east is left. The insert is a 
zoom (8x8 arcsec) around the HLX-l counterpart 
and the blue background galaxy at z = 0.03. A 
Gaussian smoothing with cr ^ 0.1 arcsec has been 
applied to the FUV image. 




Fig. 2. — Projected mass density of stars in the 
simulation at t — 2.0 Gyr after the first pericen- 
tre passage. The primary galaxy is seen edge-on. 
The scale of the color-coded map is logarithmic, 
ranging from 2.2 x 10"^ Mq pc~^ to 22 Mq pc'^. 
The frame measures 50 x 25 kpc. 



2 




Fig. 3.— Surface brightness profiles (SBPs) of the HLX-f counterpart. From left to right: F160W, F775W 
and F555W. Circles: HST data (filled circles: first approach; open circles: second approach). Dotted red 
line: stellar population component (SPC), that is simulated SBP at t = 2.5 Gyr after the first pericentre 
passage; dashed blue line: irradiated disc component (IDC), smoothed over a two-dimensional Gaussian 
profile with the same full width at half maximum (FWHM) as the point spread function (PSF); solid black 
line: total SBP from the model. Vertical dotted line: softening length. Vertical solid line: PSF FWHM. 



an annulus of width 0.08 arcsec at r > 1 arcsec 
(in the following, we refer to this procedure as 
first approach). In F555W and F775W, because 
of the high SO flux, we attempted also a different 
approach, and tried to compute the SBPs after 
subtraction of the galaxy emission (via Gaussian 
smoothing with a = 0.4 arcsec). In the following, 
we refer to this procedure as second approach. We 
do not attempt to extract a SBP from the F160W 
filter, as the contamination by the SO is too high. 

3. Model and simulations 

A^— body simulation: As in M12, the initial 
conditions for both the primary galaxy and the 
secondary galaxy in the iV— body model are gener- 
ated by using an upgraded version of the code de- 
scribed in Widrow, Pym & Dubinski (2008). Both 
the primary and the secondary galaxy have a dark 
matter (DM) halo, a stellar bulge and a stellar 
disc. The giant SO galaxy has no gas, whereas 
the secondary galaxy has an initial gas mass of 
1.38 X 10® Mq, distributed according to an ex- 
ponential disc. The total mass of the secondary 
is ~ 1/20 of the mass of the primary, classifying 
the outcome of the interaction as a minor merger. 
The orbit of the satellite is nearly parabolic, with 
initial relative speed 50 km s~^. As in M12, 
we simulate the evolution of the models with the 
A''— body/SPH tree code gasoline (Wadsley, Stadel 



& Quinn 2004). In the simulation, the particle 
mass in the primary galaxy is 2.5 x 10'* Mq and 
5 X 10^ Mq for DM and stars, respectively. The 
particle mass in the secondary galaxy is 2.5 x 10"^ 
Mq for DM and 5 x 10^ Mq for both stars and gas. 
The softening length is e = 30 pc. Fig. [2] shows 
the projected density of stars in the simulation at 
t = 2.0 Gyr after the first pericentre passage. We 
can derive synthetic fluxes from the A^— body sim- 
ulations, on the basis of each particle mass and 
age. In particular, we use the single stellar popu- 
lation (SSP) models based on the tracks of Marigo 
et al. (2008fl The tables of the SSP integrated 
magnitudes were implemented in the TIPSY visu- 
alization package for A^— body simulationtH. 

Irradiated disc component (IDC): We add 
an IDC to the synthetic profiles obtained from the 
A^— body simulation. The IDC was modelled with 
a code developed for computing the optical lumi- 
nosity of ULX binaries (see Patruno & Zampieri 
2008, for more details about the code). A stan- 
dard Shakura-Sunyaev disc is assumed and both 
the X-ray irradiation of the companion and the 
self-irradiation of the disc are accounted for. The 
model shown in Figs. [3] and |4] assumes IMBH mass 
10** Mq, bolometric luminosity L = 10"*^ erg s""'^, 
inner radius 3 (where rg is the gravita- 



^http: //stev. oapd. inaf . it/cgi-bin/cmd_2 . 3| 

^http: //www-hpcc . astro .Washington. edu/tools/tipsy /tipsy . html 
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Fig. 4.— SBPs of the HLX-1 counterpart. From left to right: F390W, F300X and F140LP. Filled circles: 
HST data, (first approach). The lines are the same as in Fig. [31 but for the three bluer filters. 



tional radius), outer radius rout — 3.5 x 10'^'^ cm, 
inclination of 45° and albedo 0.7. 

4. Discussion and conclusions 

Fig. El shows that there is good agreement 
between the observed magnitude of the HLX-1 
counterpart and the model that combines SPC 
and IDC, in the fihers F775W and F555W. In 
the F160W filter, the total magnitude of the 
model (22.87) is lower than the observed one 
(23.49±0.26). The matching in the F160W filter 
can be slightly improved (22.91 versus 23.49±0.26) 
by accounting for a realistic reddening (Ay = 
0.18). In the bluer filters (especially F300X and 
F140LP) the observed SBP shows a tail at r > 0.2 
arcsec that extends beyond the PSF limits, cannot 
be due to the SO galaxy and cannot be explained 
with the young stellar population in the simulated 
satellite (Fig.Hl). 

In summary, the comparison between the HST 
photometric data of ESO 243-49 and the results 
of our simulations confirms that a minor merger 
is a viable scenario to explain the properties of 
HLX-1 and of its counterpart. The nature of the 
extended FUV and NUV emission deserves further 
investigation, as well as the existence of tidal tails 
surrounding the counterpart of HLX-1. 
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